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Abstract: This technical report reviews the tests performed at the Best Cyclotron Systems, Inc. 
facility in regards to developing a cost effective ion source, beam line transport system, and accel¬ 
eration system capable of high current output for the IsoDAR (Isotope Decay At Rest) exper¬ 
iment. We begin by outlining the requirements for the IsoDAR experiment then provide overview 
of the Versatile Ion Source, Low Energy Beam Transport system, spiral inflector, and cyclotron. 
The experimental measurements are then discussed and the results are compared with a thorough 
set of simulation studies. Of particular importance we note that the Versatile Ion Source (VIS) 
proved to be a reliable ion source capable of generating a large amount of current. The results 
suggest that with further upgrades, the VIS could potentially be a suitable candidate for IsoDAR. 
The conclusion outlines the key results from our tests and introduces the forthcoming work this 
technical report has motivated. 
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1. Introduction 


1.1 The IsoDAR experiment 

In the presently accepted 3-neutrino oscillation model, the three mixing angles and mass 
squared splitting values are relatively well known [1]. There are, however, intriguing anomalies 
that do not fit within the three-flavor paradigm and suggest new physics beyond the standard model. 
These anomalies can be resolved with the hypothesis of a (3 -i- N) sterile neutrino model, in which 
there are three light mostly-active neutrino mass states and N massive mostly-sterile neutrino mass 
states. Among the anomalous experiments, LSND observed a 3.8 a excess in Vg events [2], Mini- 
BooNE observed a low energy Vg and Vg excess of 3.4 a and 2.8 a in a muon-neutrino flavored 
beam respectively [3], and the short baseline reactor experiments observed a Vg deficit from unity 
of 0.943±0.023 [4]. 

The IsoDAR experimental program is being developed as the first stage of the DAE5ALUS 
experiment to primarily investigate the aforementioned anomalies by addressing the sterile neutrino 
hypothesis. IsoDAR will accurately map the Vg disappearance oscillation neutrino wave within a 
1 kiloton-class liquid scintillating (ES) detector to determine the number of extra sterile neutrino 
flavors. The IsoDAR neutrino source relies on a low-energy, high-current (5 mA) Hj beam im¬ 
pinging on a ^Be target producing neutrons which flood a sleeve containing highly enriched ^Ei, 
thereby creating ^Ei (neutron capture) which decays via a well-understood j3-decay [5]. The usage 
of H 2 can mitigate the current limitations constraining the previous generation of cyclotrons and 
will enable IsoDAR to produce a scientifically useful and interesting neutrino flux. IsoDAR, paired 
with an observatory like KamLAND [6], would observe 8.2 x 10^ reconstructed inverse beta-decay 
events in five years. With this data set, IsoDAR could decisively test the N = I and N = 2 sterile 
neutrino oscillation models within the confidence intervals of the anomalous experiments, allow 
precision measurement of Vg — e^ scattering [7], and search for production and decay of exotic 
particles [8]. 

IsoDAR is split into three main components: the IsoDAR injection system, the compact cy¬ 
clotron (see Eigure |I]), and the antineutrino-production target. The latest test stand built at the Best 
Cyclotron Systems, Inc. (BCS) facility in Vancouver, Canada, was adapted to test the concept of 
the current IsoDAR injection system and is the subject of this technical report. The outcome of 
these initial tests, and this report, represent the early stages of exploration into the development of 
the IsoDAR injection system. 

1.2 IsoDAR injector requirements 

To accomplish the physics goals of the IsoDAR experiment, outlined in [9], the IsoDAR in¬ 
jector is required to address thr'ee challenges: space-charge effects, high intensity H 2 from an ion 
source, and efficient injection into the cyclotron. 

1. Space-charge effects: A limiting factor constraining the maximum current in modern day 
cyclotrons is space-charge. Self-generated electric fields due to a high charge density cause 
the beam emittance to increase, making it more difficult to efficiently inject, accelerate, and 
extract ions from a cyclotron. The generalized perveance, K, a measures of the strength of 
space-charge, is given by: 
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Figure 1. Relative size of the IsoDAR injection system and the compact cyclotron. The accelerated 
beam exits the image on the right and is transported to the antineutrino-production target approximately 
50 m downstream. 




(1.1) 


[10] with q, I, m, y, and j3 the charge, current, rest mass, and relativistic parameters of the 
particle beam, respectively and fe the space charge compensation factor. The higher the value 
of K, the stronger the space-charge effects. IsoDAR has specifically chosen to overcome 
space charge limitations. Each H 2 ion provides 2 protons after stripping thus reducing the 
necessary beam current by a factor of 2. So, accounting for the changes in mass and velocity, 
this leads to a reduction of K by a factor of 1.4 compared to protons. Space charge effects 
can be further reduced by maximizing the beam energy term), controlling space-charge 
compensation through beam line pressure (fg term), and avoiding electrostatic potentials in 
the beam line (e.g. Einzel lenses) whenever possible. 

2. Intensity of ion source: To capture 5 mA of Hj in the cyclotron, the ion source may need 
to provide up to 50 mA of Hj, depending on the combined efficiency of fhe buncher, beam 
fransporf, infiecfion, and RE capfure. While many profon ion sources readily achieve 50 mA, 
this is at the edge of what ion sources have produced so far. 

3. Injection into the cyclotron: The axial injection into the cyclotron for the more magnetically 
rigid HJ ions (compared to protons - H+), requires an unusually large spiral inflector to redi¬ 
rect the beam into the center plane of the cyclotron. The large spiral inflector is constrained 
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Figure 2. A complete schematic of the BCS test stand. The enumerated components are: [a] RF magnetron, 
[b] graded, insulated waveguide, [c] hydrogen mass flow controller, [d] quartz window, [e] plasma chamber, 
[f] high voltage side, [g] groud side and upstream ion gauge, [h] Solenoid 1 (SNl), [i] 4-Sector collimator 
(4Set), [j] water-cooled collimator (Cl), [k] first set of steering magnets on a small adapter tube (AT), [1] first 
gate valve, [m] 4-Jaw collimator (4Jaw), [n] Retarding Field Analyzer (RFA) and downstream ion gauge, [o] 
interchangeable Faraday cup (FC) and pneumatic beamstop (BS), [p] DC Current Transformer (DCCT), [q] 
Allison-type Emittance Scanner (ES), [r] second set of steering magnets, [s] Solenoid 2 (SN2), [t] 1 MeV 
test cyclotron, [u] cyclotron cryopump, [v] quadrupole focusing magnets, [w] beam line cryopump, [x] 1500 
Fs turbo pump, [y] spiral inflector (SI), [z] second gate valve. A list of commonly used abbreviations can be 
seen at the bottom of the figure as well as the approximate length of the beam line. The gray boxes represent 
the approximate locations and sizes of the important component used in the simulations (cf. Section The 
end points for Configurations I-III as described in the text are also indicated. 

by the physical space in the central region, allowable beam emittance, and maximum voltage 
across the inflector plates before electrostatic breakdown. 

The adaptations for Hj ions to the BCS test stand facility were aimed at addressing the viabil¬ 
ity of overcoming these challenges. 

2. Experimental setup 

2.1 Injector test stand layout 

The ion injector is divided into four individual sections: the Versatile Ion Source (VIS), the 
low-energy beam line transport (LEBT) system, the spiral inflector (SI), and the cyclotron. The 
experimental layout of the components can be seen in Figure @ along with a list of commonly 
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referred to abbreviations and their location along the beam line. During the runs in the summers of 
2013 and 2014, three different configurations were used to study the beam dynamics in the LEBT 
and the spiral inflector performance: 

• Configuration I: The first tests were performed with the beam line extending 3.5 m beyond 
the extraction aperture. As seen in Figure @ marked by grey boxes, the main beam shaping 
elements were a solenoid magnet (SNl) close to the source, and several sets of collimators. 
An emittance scanner (ES) and a Faraday cup (EC) were installed in two 6-way crosses 
towards the end. 

• Configuration II: The beam line was later extended to 5.5 m, with the addition of another 
solenoid (SN2) and a quadrupole doublet. Most tests using the quadrupole doublet did not 
yield any improvements for cyclotron injection and thus it was mostly unused. For the rest of 
the paper, it will be considered a simple beam pipe. The emittance scanner and the Faraday 
cup were moved to the new end of the beam line. 

• Configuration III: Finally, the cyclotron was attached to the beam line end, replacing the 
emittance scanner and Faraday cup. The rest of the FEBT is identical to Configuration II. 

Configurations I-III are labeled at the top of Figure @ and indicate the end points of the three 
beam line configurations. 

2.1.1 Versatile Ion Source (VIS) 

The VIS is a 2.45 GHz, non-resonant microwave ECR ion source that was constructed at the 
Istituto Nazionale di Fisica Nucleare, Faboratori Nazionali del Sud (INFN-FNS). It was developed 
as an evolution of the TRIPS source [11] in order to produce a robust, high current, continuous wave 
(CW) proton beam. In the VIS, ions are produced by collisions with electrons resonantly heated 
by 2.45 GHz microwaves which are supplied through an electrically insulated waveguide. The 
magnetic field necessary for elecfron cyclofron resonance is provided by fwo permanenf magnef 
solenoids. This is a convenienf configuralion as virtually no power has fo be Iransferred fo fhe 
high volfage plalform. Due fo ifs efficiency and high-currenf characferisfics, if was felf fhaf if could 
pofenfially be a good source for H^ ions, and so was shipped fo Vancouver and incorporafed info 
fhe BCS fesf sfand for fhe experimenfs conducfed by fhe collaborafive effort of fhe MIT IsoDAR 
group, BCS slaff, and fhe learn from INFN-FNS [12]. 

The VIS consisls of elemenls [a]-[g] in Figure Microwave power is provided by a Sairem 
2 kW 2.45 GHz magnelron, wifh flexible operalion in Confinuous Wave (CW) or pulsed mode. 
Diafomic hydrogen gas is injecfed info fhe plasma chamber fhrough an El-Flow F-200CV mass 
flow confroller. A 100 kV, 100 mA FuG Elekfronik GmbH power supply provides fhe high volfage 
fo fhe plalform and ion source body relalive fo fhe grounded beam line, Ihus defining fhe ion 
beam energy. Typically, due fo fhe high relalive humidify in Vancouver, fhe VIS was operafed al 
< 66 kV source polenfial, while a suppressor volfage of < —3.5 kV was used fo avoid eleclrons 
slreaming from fhe FEBT info fhe plasma chamber wifh fhe aim fo prevenl elecfrosfalic breakdown 
and discharges. 
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Figure 3. Relevant cross sections for the main physical processes in a hydrogen ion source. From [13]. 

The dominant processes leading to production and loss of in a plasma ion source are [13]: 

H 2 T ^ ^ 

H+ + e —^H+ + H + e 

! + (2.1) 

H++H 2 —^H+ + H 

H+ + e^ — >2H+ + 2e-. 

and the relevant cross sections are plotted in Figure 0. 

Atomic physics arguments show that the ratio of HJ to H+ ratio is strongly affected by electron 
density n^ and ion lifetime in the plasma T, [14]. In particular the increase of T, favors proton 
generation since it increases the probability that collide with electrons leading to their break-up 
into protons. T, can be modified in many ways, for example by modifying the magnetic field or 
fhe dimension of the plasma chamber. Since the VIS magnetic field is generated by permanent 
magnets [11], modifying the magnetic field required shifting the magnetic assembly along the 
plasma chamber. These are further investigated in Section 

The extraction system of the VIS is a typical accel-decel tetrode system with the following 
electrode voltages used during these tests: The source (plasma aperture) at +55 kV to +65 kV, the 
Puller at —2 kV to —4 kV, and the remaining two electrodes grounded. The plasma aperture had a 
diameter of 8 mm, while the apertures of the remaining electrodes were 12 mm. 

2.1.2 Low Energy Beam Transport (LEBT) 

The LEBT (elements [h]-[z]) shown in Figure is a transport system that applies dynamic 
modifications and diagnostic services to the beam. In particular, it must focus, align, and sepa¬ 
rate the ion species prior to injection into the spiral inflector. The LEBT was originally designed 
using the ray-tracing code TRACK [15], not including space charge effects. The measured beam 
parameters were later compared with WARP [16] simulations (cf. Section^. 
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Focusing of the beam is accomplished using both solenoids and a quadrupole-doublet. The 
asymmetric particle rotation combined with the radial fringe fields of the two solenoids, SNl and 
SN2 (element [h] and [s]), provide a second order focusing to the beam. The focusing length 
associated with each solenoid is proportional to the mass squared times the velocity squared, 
and inversely proportional to the charge of the ion squared times the magnetic field squared, 

The difference mmjq for profons and Hj allows fhe solenoids fo selectively focus eifher ion species 
while dispersing fhe ofher. This was fhe primary mefhod of ion species separafion wifhin fhe beam 
line. The dispersed ion species collides wifh fhe wafer-cooled collimafors (elemenfs [j] and [m]), 
or fhe vacuum lube walls. The rofalable quadrupoles (elemenl [v]), can serve lo heller malch fhe 
beam lo Ihe spiral infieclor before entering Ihe cyclolron (elemenl [w]). However, in practice Ihey 
were nol used as Ihey appeared lo have no noticeable effecl on beam transmission. 

The LEBT was aligned using a theodolite located downstream of the cyclotron, sighting with 
fiducials specifically made for various parls of Ihe beam line (e.g. axial opening of Ihe cyclolron). 
Each componenTs axis was aligned lo Ihe axis of Ihe VIS exlraclion hole lo wifhin ±0.25 mm. 
Beam skew arising from Ihe remaining misalignmenl was compensated by using dipole steering 
magnels (elemenls [k] and [r]). Two gate valves (elemenls [1] and [z]), interlocked wilh Ihe pressure 
and relevanl power supply sensors separate Ihe vacuum system into Ihree parls: Ihe VIS, Ihe EEBT, 
and Ihe cyclolron. Vacuum in Ihe upslream side of Ihe EEBT was suslained by a 1500 1/s lurbo 
pump (elemenl [x]), and in Ihe downslream EEBT and cyclolron by cryopumps (elemenls [w] and 
[u]). Hoi filamenl ion gauges were located al Ihe ion source exil (elemenl [g]), in Ihe middle of 
Ihe beam line (elemenl [n]), and in Ihe cyclolron (elemenl [1]). Experimenlal resulls perlaining to 
beam Iransporl can be found in Sections and P3| . 

2.1.3 The spiral inflector 

The axial injection of an ionized beam into a cyclolron is usually realized using an eleclroslalic 
device called a spiral inflector, which consisls of Iwo curved eleclrode deflectors (see Eigure H). The 
eleclroslalic potential belween Ihese eleclrodes is able to bend Ihe beam 90° from Ihe axial line to 
Ihe median plane of Ihe cyclolron. The helical Irajeclory of Ihe beam is determined bolh by Ihe 
shape of Ihe eleclrodes (eleclric field) and Ihe magnetic field of Ihe cyclolron as Ihe beam is benl 
into Ihe median plane. The applied vollage depends on Ihe velocity of Ihe beam as well as Ihe 
rigidity of Ihe ions. 

The spiral inflector tested al BCS was designed to mimic Ihe conditions required for Ihe Iso- 
DAR inflector, ils primary defining characteristic being Ihe 15 mm gap belween Ihe eleclrodes. 
Compared to olher similar designs [17], Ihis is ralher large because il has to lake into accounl 
Ihe larger beam size due to space-charge effecls. The spiral inflector has a copper housing lhal 
surrounds Ihe eleclrodes to minimize Ihe interaction belween Ihe eleclric fields generated by Ihe 
dees and Ihe eleclroslalic fields belween Ihe eleclrodes. To minimize beam slriking Ihe uncooled 
eleclrodes, a water-cooled, grounded, reclangular collimator shields Ihe spiral inflector enlrance. 

The preliminary design of Ihe spiral inflector called for a nominal vollage of ± 11 kV and a till 
angle of 16° for optimal beam injection al 60 keV. The lolal heighl of Ihe device was approximately 
80 mm. These parameters were found by using a MATEAB code based on an analytic Iheory for 
spiral infleclors [18,19]. The eleclrode shape was calculated using VeclorEields OPERA [20]. This 



Figure 4. CAD rendering of the spiral inflector and particle trajectories. The positive and negative elec¬ 
trostatic deflectors are shown in blue and green respectively. The cyclotron magnetic held in this image is 
directed vertically upwards. Particles enter the spiral inhector via the rectangular grounded collimator (solid 
gray) and are guided into the cyclotron mid-plane by means of the cyclotron magnetic held and the elec¬ 
trostatic potential between the electrodes. The copper housing (transparent gray) isolates the spiral inhector 
from the RF helds driving the cyclotron. The recesses on both sides of the housing (one of them visible on 
the lower left) provide space for the tips of the dees. 

design was later modified to account for some effects due to fringe fields on the electrodes. One of 
these modifications was to reduce the overall length of each electrode by 5 mm. 

First turn acceleration is achieved by the dee tips extending into recesses in the spiral inflector 
housing. The shapes of the dee tips and the spiral inflector housing have been designed to guide 
the particles from the spiral inflector exit to the acceleration region while providing the necessary 
energy gain and beam focusing. 

The final simulation (neglecting space-charge effects) of the spiral inflector was able to trans¬ 
port a 60 keV beam with a normalized emittance of 0.62 Ti-mm-mrad at a transmission efficiency 
of 100%. 

2.1.4 The test cyclotron 

The test cyclotron components include a magnetic structure, a vacuum system and an RF 
system (see Figure The magnetic structure consists of an electromagnet split in two halves, 
each enclosed in an identical steel structure. The combined steel assembly has an outer dimension 
ofl.lmW X l.lmHx 0.9 m L and opens up at the cyclotron median plane. Both halves of the 
cyclotron are mounted on rails which allow for quick and easy access to the inside of the cyclotron 
for adjustment and changes to the inflector and central region. There are five ports at which radial 
probes can be inserted to measure beam currents and profiles in the cyclotron. Four of them are 
located in the four hill regions and one in the bottom valley, this is most easily seen in Figure [T^ 
During the experiments, three radial probes were available, to be inserted at any three of these 
locations (only one is shown in Figure The inner structure is composed of permanent pole disks 
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Figure 5. Cyclotron component layout (see also Figure |T^. 


on which rest machinable components (shims and center plug) designed to replicate the IsoDAR 
cyclotron magnetic field. In order to generate the average magnetic field of 1.1 Tesla, the two large 
coils are powered by two highly stable 125 A, 80 V power supplies connected in parallel. The 
energy gain by the particles is limited to 1 MeV to avoid material activation. 

The main tank is brought down to an operating pressure of 1.33 • 10^^ mbar by a combination 
of a 605 1/min rotary vane pump and a 1500 1/s cryopump. The pressure is measured by both a 
thermocouple and a hot filament ion gauges at their respective operational range. 

Like the magnetic structure, the RF system is composed of permanent and interchangeable 
parts to achieve different cyclotron configurations. The accelerating gaps and resonator geometry 
are determined by a set of dees, center post, and center liners that are specific for each configurafion. 
Coarse funing is done manually by adjusfing shorfing plafes befween fhe sfems and fhe coaxial 
liners and fine funing is carried ouf by mofor-driven funers and couplers. The fronl-end sfage of 
fhe 20 kW RF amplifier is a wideband solid-sfafe driver direcfly connecfed fo fhe final sfage funed 
for fhe operational frequency of 49.2 MHz. The amplifier has been funed fo a maximum CW 
oufpuf power of 13 kW (17 kW equivalenf in pulsed mode). The resonators have been designed 
to reach up fo 70 kV of accelerating dee volfage. Despite our besf effort, arcing and cavify losses 
only allowed a dee volfage of < 63 kV af fhe peak power of 17 kW. The accelerating volfage 
was esfimafed using a capacifive pick-up probe. The probe was calibrafed based on measuremenfs 
of x-rays (using fhe endpoinfs of fhe bremssfrahlungs specfra) wifh a defecfor interfaced wifh a 
vacuum aperfure looking info fhe cyclofron chamber. Large flucfuafions limifed fhe value of fhese 
pick-up probe measuremenfs and if appears fhaf fhe acfual volfage in fhe cenfral-region accelerating 
gaps was even lower fhan fhe esfimafes from fhe x-ray measuremenfs. Resulfs of fhe injecfion and 
acceleration sfudies will be reported in Secfion P3| . 
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Figure 6. Left; An increased beam line pressure allowed for the beam to be easily seen. A photograph of the 
beam was taken through a quartz window situated directly below the emittance scanner at Configuration I 
(indicated in Figure 2) with 340 A supplied to SN1. The base of the emittance scanner can be seen above the 
beam. Right; An example normalized vertical emittance plot of the beam is shown. The blue dotted ellipse 
represents the 1-rms emittance and the green dotted ellipse represents the 4-rms emittance. The color bar 
represents the relative beam current and should not be interpreted as a direct measurement of total current. 
The protons and are indicated. Here, the Hj is coming into focus while the protons are highly divergent. 
What appears to be a distortion in the beam in the lower right is caused by the superposition of and 
protons. 

2.2 Beam line diagnostics 

This section describes the diagnostic instrumentation used to make measurements regarding 
beam emittances and beam currents. The data gathered using a retarding field analyzer, which 
pertains to the space-charge compensation of the beam line, will be discussed in subsequent papers, 
so is not included in this technical report. A description of the Allison-type emittance scanner, 
Faraday cup, and beamstop follows. 

2.2.1 Allison-type emittance scanner 

The kinematics of the beam are described by the particle ensemble representation in 6- 
dimensional phase space; 3-dimensions for its spatial components x, y, and z, and 3-dimensions 
for its momentum components x’, y’, and z’. The Allison-type emittance scanner (element [q]) of 
Figure is a device to measure the emittance of the beam. If {x^) is the second moment of the 
beam’s horizontal phase space distribution f{x,x')', defined as; 


2 v ^ JJx^f{x,x')dxdx' 

If f{x,x')dxdx' 

and analogous for and (xx'), fhen fhe horizonfal rms emittance can be defined as; 


( 2 . 2 ) 




[mm-mrad] 


(2.3) 


Similarly the vertical emittance can be calculated by replacing x with y and x’ with y’. The emit¬ 
tance scanner samples either the x-x’ or the y-y’ phase-space by means of two slits, electrostatic 
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Figure 7. The effect of measured beam current as a function of electron suppression voltage on the Faraday 
cup and DCCT. The current stabilizes when a negative suppression voltage of -180 V is applied. 


deflection plates, and a small Faraday cup. The spatial information comes from the position of the 
emittance scanner, actuated by a stepper motor, and the angle information from the voltage applied 
to the deflection plates located between the two slits. The current recorded in the Faraday cup for 
each pair of x and x’ gives the distribution f{x,x'). A sample emittance plot is shown in Figure ^ 
(right). 

2.2.2 Faraday cup and beam stop 

A Faraday cup is a device typically made out of copper that measures the net current incident 
on the interior face of a cup-shaped (sometimes a cone or a semi-sphere) conductor. This cup is 
connected to an ammeter or current amplifier, while an electron suppression electrode is mounted 
upstream to provide an negative electrostatic potential to turn back electrons emitted from the cup 
interior surface upon the impact of beam ions. Both cup and electron suppressor are shielded by 
a grounded aperture or housing. The tests at BCS employed a Faraday cup to accurately measure 
the total transported beam current at the end points of Configuration I and IF When the cyclotron 
was connected to the beam line, the Faraday cup was replaced by a beam stop ([o] in Figure 
that could be pneumatically moved out of the way to not obstruct the beam. This beam stop (and 
its upstream collimator) were insulated, so could also read current. Electron suppression of the 
beam stop was achieved by placing a strong permanent magnet close-by. This did not lead to full 
suppression and readings from the beam stop have a higher uncertainty (cf. Section 

2.2.2.1 Electron Suppression 

When a beam strikes the surface of a metal, electrons can be ejected from the surface. If the 
electrons are not suppressed (returned to the surface of the metal), this appears as a positive current. 
Electron suppression of the Earaday cup was accomplished by application of a negative bias to a 
dedicated electrode between the cup and the grounded entrance aperture. In order to determine the 
effect of electron suppression, a set of measurements was taken with the Earaday cup located at the 
endpoint of Configuration I and a Bergoz DC Current Transformer (DCCT) located one position 
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further upstream. As shown in Figure |7|, an electron suppression voltage of approximately -180V 
was sufficient to prevent electrons from escaping the Faraday cup. Whenever we refer to a Faraday 
cup measurement, adequate suppression can be assumed. It should be noted that an electrostatic 
electron suppression voltage was not needed for measurements on the radial probes within the 
cyclotron because the parallel magnetic field effectively suppressed fhe escaping elections. Due 
fo fhe inabilify fo suppress all fhe back-sfreaming elections artificially increasing fhe measured 
currenf in fhe DCCT, we removed fhe DCCT and unless ofherwise nofed, all subsequenf currenf 
measuremenfs were done using only fhe election suppressed Faraday cup. 

2.3 Data acquisition 

Beam currenfs were measured using custom-built bipolar current-sampling data-acquisition 
modules, connected through shielded 50 Ohm cables. The sampling rate of the input stage was 5.5 
kHz, with data averaged and stored in a buffer at a 66 Hz rate. Readout from the PLC buffer to the 
personal computer was done through Ethernet at a 5 Hz polling rate. 

Initially, upon switching to pulse-mode (required because of power limitations in the cyclotron 
RF system), it was noticed that the current amplifiers were overesfimafing fhe currenf in each pulse. 
This was found fo be due fo fhe pulse defecfion algorifhm in fhe currenf amplifiers and fhe pulse 
shape. An oscilloscope frace of several Hj pulses is shown in Figure ^ (leff). The quick rise in fhe 
pulse reflecfs fhe high initial rafe of H^ production and may be referred fo as "pre-glow" [13]. This 
can be explained by fhe production processes shown in M The production of protons is a fwo-sfep 
process involving energefic collisions wifh free elections in fhe plasma chamber, whereas fhe H^ 
is readily produced after a single process involving fhe H 2 gas. Since fhe profons producfion will 
Iherefore have a longer rise time, an initially high H^ pulse followed by an exponential decay fo 
fhe equilibrium sfate is expected. (The frace in Figure represenfs fhe Faraday cup signal when 
beam focusing is sef fo optimize for H^ fransmission fo fhe cup. When fhe beam line is sef fo focus 
profons info fhe cup, fhe pre-glow peak is missing, replaced by a slow rise-lime pulse wifh a time 
constanf of aboul 1 ms.) The pulse defecfion algorifhm in fhe currenf amplifiers would average 
over fhe inilial pre-glow and plaleau currenfs, and overesfimafe fhe acfual amounf of currenf in 
each pulse. A low-pass filter was inserfed befween fhe BNC inpuf and dafa acquisition module fo 
provide a more slable and even inpuf signal (see Figure righf). 

2.4 Uncertainty estimates 

In keeping with the goal of this paper as an examination of the concept rather than a demon¬ 
stration, we have conservatively estimated the uncertainties for the different types of measurements 
and only provide upper limits for the maximum errors. The following subsections describe these 
estimates, both statistical as well as systematic in nature, for the various types of measurements. 

2.4.1 CW beam measurement uncertainties 

The majority of the tests performed before injecting beam through the spiral inflector were con¬ 
ducted with a continuous (DC) beam. Associated with each measurement we identify an uncer¬ 
tainty associated with the measuring device as large as the least significant digit of measurement, 
as well as include a noise estimate if there was variation associated with the measurement. 
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Figure 8. Two oscilloscope traces of the beam pulse. Left; Unmodified pulse shape input into DAQ. Right; 
Capacitor inserted into DAQ board acts to average the pulse and limits the preglow pulse height. Giving a 
more accurate current measurement. 


2.4.1.1 Faraday Cup 

We assume no systematic errors on the fully electron suppressed Faraday cup as described in Sec¬ 
tion The statistical variations were typically on the order of ±2%. We conservatively report 

all beam current measurements in the Faraday cup with an uncertainty of ±5%. 

2.4.1.2 Beam stop 

The electron suppression on the beam stop was achieved by placing permanent magnets on the 
front and back of the beam stop plate, thereby creating a more or less parallel magnetic field that 
guides secondary electrons back to the plate. This cannot be considered full electron suppression 
and thus a systematic error of -10% is assumed in addition to the statistical error of ±5%. 

2.4.1.3 Beam Stop Front Plate and Spiral Inflector Aperture 

Both of these plates did not have any significant electron suppression through either a negative 
electrode or a magnetic field. They are both considered completely unsuppressed. Taking the 
curve in Figure 0 as a basis, we see that the suppressed current is ~ 60% of the unsuppressed 
one. This was measured for a copper cup with impinging HJ ions which is similar to the other 
unsuppressed surfaces. However, the Faraday cup also has a certain depth which may recover 
some of the escaping electrons. Taking this into account, we conservatively estimate the systematic 
error as -50% with an additional uncertainty of ±20% and thus report the current measured on an 
unsuppressed copper surface as factual = 0.5 • Was. ±0.2 • Imeas.- 

2.4.2 Pulsed beam measurement uncertainties 

As discussed in Section P3| , the pulsed beam exhibited a pre-glow effect. A capacitor was 
added, effectively working as a low-pass filter mostly evening out the peak at the beginning of 
the pulse. Comparison of measured values with oscilloscope traces lead to estimating the system¬ 
atic error for any pulsed beam as ±5% from this effect, in addition to the previously discussed 
uncertainties for the various devices. 

2.4.2.1 Paddle and Radial Probes 

Inside the cyclotron a strong magnetic field on the order of 1 T was present which leads to full 
electron suppression during the current measurements. We thus report the measured currents with a 
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base statistical error of ±5%. Additional systematic errors from probes picking up RF signals from 
the dees during acceleration tests are taken into account by subtracting the background obtained by 
inserting the beam-stop and measuring the signal without ion beam in the cyclotron. 

2.4.3 Emittance measurement uncertainties 

The uncertainties in the reported measurements come from the quadrature combination of un¬ 
certainties in the measuring devices and observed electronic noise. Further, we add an uncertainty 
of ±5% due to the electron suppression described in Section [2.2.2.11 There is an additional large 
systematic error from emittance scans not covering the full phase space area of the beam. This 
systematic error can be determined by simulations and is included in the values reported in the text. 
However, the graphs in Figures |T^ and 0 do not include the systematic modification. Instead, the 
simulations were restricted to the same limits, showing good agreement. 

2.4.4 Dee voltage uncertainties 

Due to the highly unreliable pick-up probe for measurements of the dee voltage and large 
fluctuations in the dee voltage itself, we do not estimate an uncertainty for the dee voltage. Instead 
we report the values as they were read, but suggest that the actual values are more closely related 
to the best fit obtained from the OPERA compaiisons to measured transmission efficiencies (cf. 
Section R3|) . 

3. Experimental measurements 
3.1 Overview 

In March of 2013 the VIS source, high-voltage platform and supply, microwave generator, 
a large-bore solenoid plus controls, power supplies and miscellaneous supporting equipment were 
shipped from Catania, Italy to the BCS facility. The shipment was met by a team from the INFN- 
LNS, who unpacked the equipment and interfaced it with the heavy-duty rails supporting the beam 
line equipment and the cyclotron. During the summer of 2013 initial tests were performed by our 
MIT group, staff of BCS and the Catania group. The first results were reported in [21]. Based 
on these results, the months between the summers of 2013 and 2014 were used to perform im¬ 
provements to the data-acquisition system, beam diagnostic devices, cyclotron RF system, and (as 
described below) to the VIS. The results summarized in this Section were mainly obtained in the 
2014 run. 

To begin addressing the IsoDAR injector requirements enumerated in Section 113, the tests 
primaiily sought to: 

1. characterize and optimize the VIS for Hj operation (Section ^3 ), 

2. separate the ion species in the transport line (Section P3[ ), 

3. determine transport efficiencies of the beam line and the spiral inflector (Section P3 ), 

4. capture beam into the cyclotron (Section P3[ ), 


- 15 - 





Figure 9. The plasma chamber designs. Left; The standard chamber; 10 x 5 cm (length x radius); Right; 
The small chamberl; 10 x 2.5 cm. 

5. verify accuracy of beam line simulations (Section ^3). 

3.2 Characterizing the VIS 

The results obtained in 2013 showed that the VIS could provide good currents. However, 
as the source design was optimized for proton currents, as expected suffered in comparison 
to proton currents. As one raised the microwave power, a stable discharge first occurred at about 
400 watts. At this point the total extracted beam current was about 30 mA, of which 8 to 10 mA 
was HJ , the remainder being protons. As the microwave power was raised, the proton current 
increased, up to a maximum of about 40 mA at 1400 W, however the current remained constant 
at 8 — 9 mA. The higher power increased the ion production, but also increased the breakup of Hj 
into protons. 

Due to the fact that the VIS is using a permanent magnet structure which was optimized for 
proton production, the increase of fraction was achieved solely by reducing the ion confinement 
time through reduction of the plasma chamber diameter. While the original cylindrical plasma 
chamber had a diameter of 10 cm and a length of 10 cm, the new one has a diameter of 5 cm 
still with a length of 10 cm. This new chamber, and a revised microwave coupling section, were 
brought to Vancouver midway through the summer of 2014. In the companion paper [12] these 
modifications are described, and results obtained on the BCS test stand are reported. To summarize; 
As expected, the fraction increased to almost 50% at low microwave powers (about 150 watts 
now), and the best current achieved now improved to about 12.2 mA. Small adjustments in the 
location of the permanent magnet ring were performed to optimize the Hj current value. Additional 
adjustments to the plasma chamber shape might improve HJ currents further. 

All beam line measurements reported in this paper, as well as all the cyclotron injection studies 
were performed with the original, standard plasma chamber. 

3.3 Ion species separation 

As the BCS test stand was initially designed to transport ions (for the development of a LEBT 
lines for isotope-producing cyclotrons), there were no provisions for magnetically analyzing the 
beam produced by the VIS. In particular, to measure the current of an ion species, a Faraday cup 
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Figure 10. Configuration I, measured and simulated beam curTents as a function of SN1 current. The beam 
current is measured in an electron suppressed Faraday cup located behind a 56 mm aperture at the end point 
of Configuration 1. Extraction voltage was maintained at 60/-3.2 kV, RF power at 400 W, and H 2 gas mass 
flow controller at 46%. 

and emittance scanner were used. Since the VIS was mounted straight along the axis of the rails, 
the protons, H^, and perhaps a small fraction of all propagate along the same axis from the 
source towards the cyclotron. Typically, in order to separate the ion species, a dipole magnet would 
be used but due to space constraints an alternative method was required. Simulations indicated that 
a solenoid with a longitudinal magnetic field strength of roughly 2.7 kG, placed 50 cm downstream 
from the source extraction aperture, could effectively be used to separate the various ions. As seen 
in Figure [T^, the measured current on the Faraday cup at the Configuration I end point agreed well 
with the beam line simulation; the relative current at both the proton focus (SN 1: 240 A) and at the 
HJ focus (SNl: 340 A), was in good agreement. 

3.3.1 Beam shape studies, formation of a hollow beam 

Figure 0 shows the beam striking the vertically separated water-cooled 4-Jaw collimator 
(4Jaw) (element [n] of Figure An intensity profile through the diameter of the beam is super¬ 
imposed onto the image. The left plot shows the transverse beam profile with 240 A supplied to 
SNl. At this setting, the protons are directly focused onto the 4Jaw, 2.30 m away from extraction 
hole of the VIS (see Figure ^ for the simulated beam envelop near this setting). At the center of 
the concentric rings is the focused proton beam. The high charge density of the protons repel the 
HJ, creating a circular halo of - The right plot shows the scenario for which SNl is supplied 
with 320 A and the is near its focus (see simulation in Figure The strong magnetic field 
on SNl causes the protons to be over-focused. This creates a region upstream of the 4Jaw with an 
extremely high positive charge density which repels the Hj. An depletion zone is then created 
at the center of the beam. The formation of this “hollow beam" was also observed in the beam 
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Figure 11. Transverse beam profiles observed on the 4-jaw slits (4Jaw). The beam enters the image from 
the right and strikes the face of two horizontally separated water-cooled copper plates. Left; The focused 
proton beam is seen as the central illuminated circle while the ions form a concentric outer ring. Right; 
The bright inner ring is the being brought into focus. By focusing the Hj, the protons are over focused 
and blast out a hole in the center of the beam, seen as the central intensity deficit. Note the horizontal splice 
in the image corresponds to the small separation between the water-cooled plates. 


line simulations and can be compared directly with the results described in Section 

3.4 Beam transport 

3.4.1 Intensity fluctuations 

Beam line pressure fluctuations were found to affect the beam current measurements. A video 
of the beam, captured through the quartz observation window, indicated that the intensity of the 
beam fluctuated in sync with the 1.25 Hz compression cycle of the helium cryopump. A Lomb- 
Scargle periodogram analysis [22] of beam line pressure data over a period of 5 minutes indicated 
that there were two significant oscillation frequencies in pressure; one at 1.25 Hz and another at 
5.05 Hz, see Figure [T^ The 1.25 Hz frequency was found to corresponded to the compression cycle 
mentioned above, indicating that the cryopump sorbent was saturated by the pumped out hydrogen 
gas, while the origin of the 5.05 Hz frequency was undetermined. To avoid downtime caused by the 
need to regenerate the cryopump, the beam line cryopump was replaced with a turbo pump capable 
of achieving the same beam line pressure. This modification lead to more stable beam transport 
across the LEBT. 

3.4.2 Emittance measurements 

Two sets of emittance measurements are presented. The first will describe the vertical emit¬ 
tance measurements made at the end of Configuration I for a transported H^ beam. The second set 
include both horizontal and vertical measurements taken at the end of Configuration II, and are of 
particular interest because of their significant aberrations. 

For all the subsequent emittance measurements, the proton contamination has been subtracted 
via software in order to obtain a more accurate representation of the H^ emittance. 
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Figure 12. Lomb-Scargle periodogram of the beam line pressure. Two significant periodicities in the data 
are seen at 1.25 Hz and 5.1 Hz. The 1.25 Hz periodicity was identified as the compression cycle of the 
cryopump. The 5.1 Hz periodicity was not identified, but was assumed to also be attributed to the cryopump. 

3.4.2.1 Measurements of Configuration I 

The Allison-type Emittance Scanner (ES) was placed in the beam line at the end point of Con¬ 
figuration I in order to characterize the beam shape before final focusing and injection into the 
cyclotron. Eigure |T^ shows the phase spaces of the beam for SNl = 310 A (top left) to SNl = 350 
A (bottom). The ion source high voltage potential was reduced to 55 kV (55 keV ion beam kinetic 
energy) in order to measure the phase spaces of the beam past its focal point. At 55 keV, the 
beam is focused onto the ES entrance aperture between SNl = 330 and 340 A. The observed mini¬ 
mum size of the beam was at 330 A, where the 1-rms beam diameter was measured to be 15.2±0.5 
mm. The 4-rms emittance measurements vary from 0.81 to 1.07 Ti-mm-mrad (0.94— 1.07 includ¬ 
ing the systematic error discussed in Section The rather strong aberrations in phase space and 
the resulting large emittances and large variation in emittance is believed to be the interplay of the 
following effects: 

• By changing the focusing strength of SNl, less and less beam is scraped on upstream colli¬ 
mators (this is directly seen in Eigure [T^ as an increase in intensity). 

• Protons and are focused differently due to their respective magnetic rigidity. 

• Eilling a solenoid magnet too much can introduce spherical aberrations into the beam phase 
space. 

• The hollow beam effect. As the protons are focused further and further upstream, they form 
a tighter and tighter spot and the space charge of the proton beam increases significantly. 
This changes the dynamics of the ions of the beam. This can be observed qualitatively 
in Eigure |T^ and is in agreement with the observation of a hollow beam discussed in Section 
|3.3.1| . In the figure, fhis is mosf evident in the top left image where there appears to be an 
HJ deficit at yss —5 mm and two Hj lobes on either side (one at y« —15 mm and the larger 
lobe at 10 mm). 
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Figure 13. Configuration I: vertical phase space measurements at several different SNl currents. The source 
potential and electron suppression were held at 55.0±0.1 kV and —3.0±0.1 kV, respectively. The RFpower 
was maintained at 400 W. The z-axis represents the measured current per pixel and was not normalized to 
the total beam current. 


In Section |4.2.1| the phase space measurements at this location will be compared to simulation 
furthering our understanding of the effects. 
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Figure 14. Configuration II; Emittance and beam diameter as a function of the current on SN2. It should 
be noted that, as discussed in Section |2.4t , there is an additional systematic error due to the fact that many of 
the measured phase spaces used to calculate the emittances did not cover the full beam area thus leading to 


artificially reduced emittance values (most notably at the beginning and end of the scan). The graph shows 


the values without the systematic errors and the simulated phase spaces, truncated with the same limits, yield 
good agreement in emittances and diameters. The systematically corrected values are reported in the text. 

One particular result of these effects is that the proton beam (not shown) does not exhibit 
similarly strong aberrations. It can therefore be concluded that the emittance of the beam 
could, in theory, be kept much lower, if one used a dipole magnet for separation of the ion species 
thereby avoiding the over-focusing of the protons. 

3.4.2.2 Measurements of Configuration II 

An ISO-160 4-way cross containing the emittance scanner (ES) and a Faraday cup was introduced 
at the end point of Configuration II (cf. Figure @) in order to take a series of transverse emittance 
measurements to ensure adequate beam quality to inject into the cyclotron. For these measure¬ 
ments, the RF power supplied to the plasma chamber was held constant at 300 W and the extraction 
potential and electrode bias were held at 62.7 kV and -3.0 kV respectively. SNl was maintained 
at 353 A in order to transport the maximum amount of to the ES. The downstream ion gauge 
measured a beam line pressure of 7.7 ± 0.1 x 10^^ mbar. The culmination of these settings yielded 
a total H 2 beam current of 5.5±0.1 mA. The current on SN2 was incrementally increased from 215 
A to 250 A. The resulting emittance measurements are summarized and Figure [T^ The horizontal 
and vertical emittance plots from SN2 = 215 to 230 A are shown in Figure |T^ and SN2 = 240 to 
250 A in Figure |T^. 

The left and right columns of Figures |T^ and |T^ show the horizontal and vertical phase space 
plots respectively, ranging from SN2 = 215 A to 250 A. Similar to the emittance plots at the end 
point of Configuration II, the thin streak seen in each graph represents the remaining proton con¬ 
tamination. For the emittance measurements recorded here, the proton contribution was subtracted 


-21 - 















from the each image to give a more accurate representation of the Hj beam emittance. As the 
current is increased on SN2 the beam diameter decreases. At 250 A, the measured 2-rms beam 
diameters along the x and y axis are « 8 mm and « 9 mm, respectively. The measured normalized 
4-rms emittance remained between 0.96 and 1.23 Ti-mm-mrad, with the 4-rms ellipses containing 
between 92.6% — 94.1% of the beam (cf. Figure [T^. Including the systematic error from cutting 
away part of the beam, these values increase to 1.15 — 1.23 Ti-mm-mrad. With the systematic error, 
the variation goes down and emittance is largely conserved for the different solenoid currents. The 
beam quality measurements shown here met the requirements for clean transport through the spiral 
inflector and provided insight into the transport of the beam through the LEBT. A discussion of the 
origin of the aberrations is presented in Section |4.2.2| . 
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Figure 15. Configuration II; a series of emittance measurements at various currents supplied to SN2. The 
horizontal emittance plots are shown on the left; the vertical emittance plots are shown on the right. The SN1 
current Ei, SN2 current I 52 , and percentage of beam contained in the 4-rms contour (green dashed contour) 
are shown in the top right of each image. The 1-rms beam diameter and 4-rms emittance measurements are 
shown at the bottom right of each plot. SN2 currents 225 A and 235 A. It should be noted that the emittances 
are without the systematic error discussed in Section [2.4| . Confer to the text for the corrected emittances. 
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Figure 16. Continuation of Figure |T^ SN2 currents from 240 A to 250 A. 
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3.5 Cyclotron injection studies 

The injection into the cyclotron was tested in three steps: 

1. Injection into an axial Faraday cup at the center of the cyclotron 

2. Injection through the spiral inflector, with measurement of the beam current on a ‘paddle’ 
radial probe. For this measurement, the central region was mounted inside the cyclotron 
rotated by 90°. 

3. Injection through the spiral inflector and acceleration. The currents are measured using three 
different radial probes. 

The three steps will now be discussed in the following subsections. 

3.5.1 Injection into an axial Faraday cup inside the cyclotron 

The first injection test was simply to replace the diagnostic box at the end of the beam line 
with the cyclotron as seen in Figure 0. Instead of the central region containing the spiral inflector, 
a Faraday cup was mounted axially inside the cyclotron. The reason for this setup was to exactly 
determine the amount of beam current entering the cyclotron through the spiral inflector aperture. 
To this end, the Faraday cup was outfitted with a similar aperture at the same z position. Electron 
suppression on the Faraday cup was provided by a suppression electrode at -180 V. The cyclotron 
magnet was engaged to 100 A 45% of the nominal current for acceleration. This improved 
beam transport into the Faraday cup significantly. The result of this first step is summarized in 
Table |T]. The transport efficiency from the beam stop (end of Configuration 1 in Figure @) to the 
Faraday aperture is close to 100% (with the current on the aperture being the main source of error 
due to missing electron suppression) and the current measured inside the cyclotron was 7.25 ±0.36 
mA leading to a final transport efficiency from beam stop into the future spiral inflector of 89%. 
The purity in the Faraday cup was determined by simulation to be > 98%. 


Table 1. First injection test with axial Faraday cup inside the cyclotron. Parameters and results. 


Parameter 

Value 

Species 

H+ 

Beam Energy 

62.7 keV 

Solenoid 1 

356 A 

Solenoid 2 

250 A 

Cyclotron Magnet 

100 A 

Beam Stop Current 

8.2l?j mA 

Earaday Cup Collimator Current 

0.5 ± 0.2 mA 

Axial Earaday Cup Cun'ent 

7.25 ±0.36 mA 
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Figure 17. A photograph of the inside of the test cyclotron. The paddle probe reaches in from the bottom 
(currently retracted to see the burn marks) to measure beam current transported through the spiral inflector 
(inside copper housing) but not yet accelerated. 

3.5.2 Injection through the spiral inflector onto the paddle probe 

The axial Faraday cup inside the cyclotron was replaced by the central region piece with the 
spiral inflector, but rotated by 90° to accommodate a paddle probe reaching into the dee recess from 
below (cf. Figure |T7|. This way, beam exiting the spiral inflector only traveled about 45° before 
hitting the paddle probe, providing a good measure of spiral inflector transmission. Secondary 
electrons from the paddle probe were suppressed by the strong magnetic field of the cyclotron. 
The currents on the spiral inflector entrance aperture were unsuppressed and as per Section 
we are reporting the readings as factual = 0.5 • Was. ±0.2 • Was.- The result of this second step 
is summarized in Table 0. For this test, the beam was pulsed and thus the additional error from 
the preglow effect combined with the low pass filter of the DAQ system is applied as discussed 
in Section Combining the measured values in Table @ with the discussed errors, we obtain 
a transmission of 93.5+J 5 g% through the spiral inflector without subsequent acceleration. This 
compares well to measurements in the next section and simulations in Section 


Table 2. Injection test with the spiral inflector in place but rotated by 90° to accommodate a paddle probe. 
Parameters and results. 


Parameter 

Value 

Species 

H+ 

Beam Energy 

60.0 keV 

Solenoids 1/2 

340 A/240 A 

Cyclotron Magnet 

223 A 

Spiral Inflector Upper/Lower Electrode 

-10kV/+10kV 

Beam Stop Current 

7.5 ±0.8 mA 

Spiral Inflector Aperture Current 

1.3 ±0.5 mA 

Paddle Probe Current 

5.8 ±0.4 mA 
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Figure 18. A photograph of the inside of the test cyclotron, the dees are labeled "D" and the radial probes 
are enumerated from (1) to (3). The beam is entering the cyclotron axially (towards the reader) is bent into 
the mid-plane and then follows a trajectory similar to the red spiral. 


3.5.3 Injection through the spiral inflector with acceleration 

The final setup for injection and acceleration tests is shown in Figure Beam enters the cy¬ 
clotron axially, is bent into the mid-plane and follows a counter-clockwise spiral trajectory through 
the cyclotron. Two dees with two accelerating gaps can be seen in the image (denoted by the letter 
’D’). During the tests, a function generator was used to pulse the ion source magnetron and the 
cyclotron RF amplifier with a 10 ms ON / 10 ms OFF pattern seen in Figure ||. Three radial probes 
where available (labeled 1 — 3 in the figure). These probes could be moved via stepper motors. 
The position accuracy could be determined to be better than ±1 mm. Readout from the probes was 
done using the current readout channels and methods discussed in Section with low pass filters 
attached to filter out the preglow peak. Because the RF amplifier could not provide the necessary 
power to obtain the full dee voltage of 70 kV, the experiments were conducted at lower dee voltages 
of ss 50kV which consequently led to lower injection efficiencies as the design was optimized for 
70 kV. A comparison with OPERA simulations of this reduced dee voltage is presented in Section 
P~3| . In addition, the RF was unstable and discharges in the resonators occurred frequently. Large 
uncertainties are therefore attached to the presented numbers. In the following, two sets of mea¬ 
surements are presented using the radial probes: Set 1 and Set 2. In these sets, the beam current 
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Figure 19. Probe current as a function of turn number for two dee voltages. Turn 0 starts at the exit of the 
spiral inflector. The parameters of the measurements are listed in Table 0. 

on each probe was measured, then the probe was retracted until no beam was seen on the probe 
anymore, indicating it was now just outside of the beam radius, thus allowing the beam to circulate 
to the next probe. This was done with all three probes sequentially until the maximum range of 
each probe was reached. The important parameters of Sets 1 and 2 are listed in Table ||. In both 
sets, the dee voltage was estimated using a capacitive pick-up probe. The values fluctuated between 
50 and 60 kV. Because of the low reliability of the probe, it is believed that the actual average dee 
voltage was on the lower end of this interval, which fits well with the simulations presented in Sec¬ 
tion P3t The measured currents as a function of fractional turn number can be seen in Figure [T^ 
the magnetic field of the cyclotron and the spiral inflector voltages were optimized individually, but 
came out to be the same for both sets. It can clearly be seen that the acceptance of the cyclotron 


Table 3. Parameters for inflection and acceleration studies. 


Parameter 

Setl 

Set 2 

Species 

H 2 + 

H+ 

Initial Beam Energy 

62.7 keV 

62.7 keV 

Solenoids 1/2 

350 A / 233 A 

350 A / 233 A 

Cyclotron Magnet 

218.6 A 

218.0 A 

Spiral Inflector Upper / Lower Electrode 

-10.0kV/-rl0.15kV 

-10.0kV/-rl0.15kV 

Beam Stop Current 

6.8±0.8 mA 

6.9T0.8 mA 

Spiral Inflector Aperture Current 

1.0±0.4mA 

l.l±0.4mA 

Approximate dee voltage 

Ri47kV 

Ri50kV 

Captured Beam 

Ri 0.4% 

« 1% 
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Figure 20. Gaussian fits to the measurement of cuiTents on the 5-finger probe. Radial FWHM of the center 
finger (Finger 3) is « 2.85 mm, vertical FWHM of the beam is « 8.2 mm. The tails on the left side of each 
data set are due to ions not fully accelerated thus having a smaller radius. 


drops with dee voltage deviating more from the design voltage. The beam radii and estimated beam 
energies are listed in Table ^ With a five-finger probe in place of Radial Probe 1, it was possible 
to show that the beam was centered vertically in the cyclotron mid-plane. With the five-finger 
probe, the beam size at the position of Radial Probe 1 during the first turn could be estimated to 
be 2.85 ± 0.5 mm FWHM radially and 8.2 ± 0.7 mm FWHM vertically (see Figure The final 
result of the injection and acceleration studies was that 66 /f A of Hj could be successfully injected 
and accelerated to 3 3/8 turns which corresponds to « 410 keV/amu and an injection efficiency of 
1% which was expected for the dee voltage of « 50 kV (cf. Section EH)- 


Table 4. Probe 1 radii, currents, and approximate beam energy for the two sets of inflection tests. The 
approximate energy gain per gap (2 per dee) was estimated to be 47 keV for Set 1 and 50 keV for Set 2. 



Approx. Radius 

Current 

Approx. Energy 

Average Mag. field 

Setl 

74 mm 

3314 AtA 

157 keV 

1.1 T 


105 mm 

109 AtA 

345 keV 

1.1 T 


135 mm 

33 ixA 

533 keV 

1.1 T 

Set 2 

75 mm 

3564 nA 

163 keV 

1.1 T 


105 mm 

231 nA 

363 keV 

1.1 T 


135 mm 

86 /lA 

563 keV 

1.1 T 


165 mm 

66 /lA 

763 keV 

1.1 T 
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4. Simulations 


One very important aspect of the design and commissioning of a particle accelerator system 
is beam transport simulations and their comparison with measurements. For the test runs at BCS in 
Vancouver, the particle-in-cell (PIC) code WARP [16] was used to simulate the transport of two ion 
species (protons and H^) simultaneously through a model of the beam line used during the tests 
(cf. Figure @). The different beam line elements and how they were included in the simulations are 
described in Section 14. 1 .21 . 

In high intensity, low energy beams, it is important to include space charge effects as well as 
make some assumptions about space charge compensation. The way this is done in the WARP 
simulations is described in more detail in Section |4.1.41 One important parameter in the compen¬ 
sation estimation is the pressure along the beam line, because it determines the density of neutral 
gas molecules available for secondary ion and electron production through residual gas ionization 
and charge-exchange processes. Simulations using the code MOLFLOW were performed in order 
to obtain pressure distributions along the beam line for several important cases and are presented 
in Section 14.1.31 

4.1 Introduction to the WARP LEBT simulations 

4.1.1 Basic simulation parameters 

WARP includes a multitude of different field-solvers which are discussed in more detail else¬ 
where [16,23]. For the presented problem, the XY-slice solver was chosen. Here, the beam is 
transported thr'ough a series of transversal slices along the z-axis (typically 0.5 mm step size) and 
the transversal self fields are calculated on a 2D grid (typically 512x512) at each step. This is a 
suitable approach for a DC ion beam in which the longitudinal self fields are largely negligible. 
The XY slice solver is robusf and fast, which makes it an ideal candidate to map a large parame¬ 
ter space. The number of macro-particles used in these simulations was 60,000 for each species 
(individual species currents are set by charge/macro-particle). 

4.1.2 Beam line components in the simulations 

The different configurations (I-III) used during the measurements were discussed in Section 
P7T| and are shown in Figure The important beam shaping elements were the two solenoids and 
the various apertures and collimators. 

The fields of solenoids SNl and SN2 were calculated in POISSON Superfish [24] and saved 
as axially symmetric 2D field maps for WARP import. Comparison of the calculated field maps 
with measurements using a Gauss-meter yielded excellent qualitative agreement. A scaling factor 
was introduced to match the peaks fields quantitatively. As can be seen in Figure |T^ this lead to 
good agreement of the proton and focal points. 

The numerous apertures and collimators were introduced in WARP as beam scrapers and care 
was taken that variable apertures (i.e. 4-Jaw slits) were set correctly. 

For simulation purposes, the LEBT in Configuration III is essentially the same as in Config¬ 
uration II with particle distributions saved at slightly different z locations. The simulation of the 
cyclotron injection was done with Vector Fields OPERA using the particle distributions obtained 
from the Configuration II WARP simulations and are treated in Section 


- 29 - 









Figure 21. Pressure distributions along the beam line in Configuration 1 calculated with MOLFLOW. It 
can be seen that the major contribution comes from the ion source gas supply with an additional pressure 
increase downstream from beam striking elements of the beam line and thereby heating and sputtering them. 
Although the pressure distributions for proton and H J foci are slightly different, they result in identical beam 
dynamics during the simulations. Also noteworthy is the effect of collimator (Cl) and adapter tube (AT), 
both with small diameter, on the pressure distribution. They mainly constrict the flow and for the different 
scenarios, different parts of the beam line contribute more or less to the overall pressure (e.g. beam induced 
out-gassing on the downstream side). 


4.1.3 Beam line pressure profile simulations 

In order to estimate space charge compensation in the WARP LEBT simulations, axial pressure 
distributions along the beam line were calculated for several different cases using the MOLFLOW 
[25] simulation package. These simulations were matched to recorded upstream and downstream 
ion gauge readings by varying desorption coefficients 

1. along the beam line (outgassing and leaks), 

2. of the ion source aperture (main gas inflow from the source). 

Three simplified models of the beam line, corresponding to the three configurations listed in Sec¬ 
tion |4.1.2| were created in the CAD software AutoDesk Inventor [26] and then imported into 
MOLFLOW. Configuration I will be used as an example in this section. 

MOLFLOW is a simulation package developed by CERN that uses a Monte Carlo algorithm 
to calculate characteristics of a system under High and Ultra-High Vacuum (HV and UHV) con¬ 
ditions. This pressure regime is called molecular flow regime where collisions between particles 
are unlikely and where collisions with the walls dominate. The molecular flow regime is reached 
when the mean free path of the gas molecules becomes longer than the vacuum vessel dimensions. 
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At 1 • 10^^ Torr, the mean free path is on the order of a few meters and so the condition is fulfilled. 
This makes Monte Carlo simulations well suited to studying a beam line under HV/UHV condi¬ 
tions. MOLFLOW takes an input geometry and generates particles according to the user defined 
desorpfion coefficienfs on fhe facefs. Pumps are defined by fhe user by absorpfion coefficienfs on 
fhe corresponding facefs. Parficles are fhen moved fhrough fhe sysfem by ray-fracing, having a 
chance of reflection or absorpfion every time fhey hif a surface. The pressure on a selecfed facef 
can be calculafed from fhe number of hifs per unif time and fhe area of fhe facef. 

The algorifhm of obfaining on-axis pressure disfribufion for fhe LEBT is as follows (confer 
also fo Figure ^ for fhe differenf sfages of fhis for Configuralion I): Inifially, fhe simulafion is run 
wifhouf any gas flow and beam from fhe ion source, so fhaf desorpfion coefficienfs corresponding 
fo leaks and room-femperafure oufgassing could be esfablished (fhe baseline). Nexf, fhe desorption 
from fhe ion source aperfure face! was increased until fhe upsfream ion gauge reading corresponded 
fo fhe measured value (fuming on fhe source gas buf nof fhe plasma). As expecfed, fhe downsfream 
ion gauge value immediately wenf up fo a value very close fo fhe reading during fhe same sfage in 
fhe experimenf. To evaluate fhe sysfem wifh beam, preliminary LEBT simulafions using a consfanf 
pressure were used fo esfablish facefs fhaf were sfruck by fhe beam (cf. Figures ^ and and 
desorpfion coefficienfs were increased again until fhe fwo ion gauges values mafched fhe measured 
ones. If should be noted, fhaf fhe gas flow from fhe source was reduced slighfly, fo accounf for fhe 
pumping effecf of fhe plasma (wifh fhe same gas flow info fhe source, fhe gas leaking fhrough fhe 
exfracfion hole will be less if plasma is ignifed). 

For fhis beam line, if fumed ouf fhaf fhe variation in on-axis pressure was small enough fo be 
secondary fo ofher effecfs (inferplay befween fhe fwo ion species, change in beam radius), buf for 
longer beam lines wifh overall lower pressures and larger variations, fhis will be an imporfanf fool 
fo esfimafe space charge compensafion. 

4.1.4 Space charge compensation in the WARP code 

Space charge compensafion arises from fhe inferacfion of beam ions wifh fhe residual gas 
molecules. There are several processes fhaf can occur, buf fhe fwo dominanf ones are charge- 
exchange and ionizafion. Togefher, fhey lead fo low energy secondary ions and elecfrons inside fhe 
beam envelope. Because of fhe posifive beam pofenfial, elecfrons are affracfed and posifive ions are 
repelled. This can lead fo a significanf lowering of fhe beam pofenfial. In fhe LEBT simulations, 
space charge compensafion is calculafed dynamically af each lime slep. The model is based on fhe 
work of Gabovich el al. [27-31] and was recenlly reviewed and updated [32,33]. In fhe model, 
fhe energy balance of fhe secondary elecfrons is used fo calculale a sleady-slale value of fhe space 
charge compensation factor fg. If depends on fhe neulral gas density (beam line pressure), fhe 
beam currenfs, fhe beam energy, and fhe beam size. If also relies on knowing fhe cross seclions for 
folal secondary ion and secondary eleclron production fhrough fhe aforemenlioned processes. This 
model is a ‘besl case’ approximalion, where fhe beam is uniform and round and all species have fhe 
same radius. No collecfive effecls, plasma oscillalions, and non-linearilies are laken info accounf. 
In many cases, however, if works remarkably well. For fhe presented applicalion, fhe difference 
in beam size had fo be laken info accounf phenomenologically. If we consider for a momenl 
fhe fwo species as independenf sysfems, we can immediately see fhaf if one of fhem is smaller in 
diameter and has a higher beam pofenfial, if acfs like a charged eleclrode for fhe second species and 
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Figure 22. Example of space charge compensation along the beam line for Configuration II. The smooth 
dips after the 4-sector collimator (4-Set) and at the electrostatic emittance scanner (ES) are from the beam 
coming to a focus after SNl and SN2. The sharp changes are from a phenomenological approach to multiple 
species with different radii and are explained in the text. 



VIS SNl 4Set Cl AT 4Jaw C+BS SN2 ESFC 


will collect electrons that would normally contribute to the second species’ compensation. In the 
presented simulations, a threshold ^ was introduced so that if both of the following two conditions 
are satisfied simultaneously: 


r\ < ^■r 2 

h > h 


(4.1) 


the space charge compensation factor fg is reduced by a factor 7 ]: 


fe.new — t] ' fe.old (■4-2) 

Here, ru are the 4-rms beam radii, I 12 the current densities, and fe.new the new compensation 
factor of the larger species. In this approach, ^ and T] are treated as free parameters. An example 
of the resulting compensation factors can be seen in Figure ^ It becomes immediately clear that 
this is a crude approximation leading to (nonphysical) sharp changes in fg. A smoother model or a 
more involved analytical solution taking into account different beam radii is desirable but not yet 
developed. For the BCS LEBT simulations, the values of ^ = 0.667 and T] = 0.8 were found to 
lead to good agreement with experiment. 

In addition, care has to be taken such that electrostatic elements in the beam line are taken 
into account. These can negate space charge compensation completely as they may collect the 
compensating electrons. 
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Figure 23. Raw initial particle distributions used in the WARP simulations, obtained with KOBRA-INP. 
The red markers represent protons, the blue markers represent Hj. Left; Initial beam cross-section. Right: 
initial x-x’ phase space. There are initial asymmetries, which can be attributed to numerical effects in the 
3D simulations. 


4.1.5 Initial particle distribution at the VIS source 

The initial particle distribution was obtained by the Catania group using the self-consistent 
3D ion source extraction simulation software KOBRA-INP [34]. The beam profile and x-x’ phase 
space 14 cm after the extraction aperture of the source are shown in Figure ^ This is the starting 
point of the WARP simulations. An extraction voltage of 60 kV was used to obtain this initial 
distribution and the longitudinal velocity component (v^) of each particle is scaled appropriately in 
the WARP simulation if a different beam energy is desired (i.e. if the experiment the simulation 
should be compared to was at a different extraction voltage). The initial 2-rms diameters and 4-rms 
normalized emittances are listed in Table 0. 

As a benchmark and to make sure the initial asymmetries were not causing the effects at¬ 
tributed to the beam line elements, a Gaussian beam with the same Twiss parameters was gener¬ 
ated and occasionally used instead of the KOBRA-INP initial distribution. The results exhibited 
the same aberrations and effects. 


Table 5. The initial diameters and emittances for protons and obtained by KOBRA-INP extraction 
simulations and used in the WARP beam transport simulations. The 4-rms emittances include « 85% of the 
beam. 



protons 

H2+ 

x-diameter (2-rms) 

15.2 mm 

15.0 mm 

y-diameter (2-rms) 

15.4 mm 

15.0 mm 

x-x’-emittance (4-rms, normalized) 

0.54 TT-mm-mrad 

0.36 TT-mm-mrad 

y-y’-emittance (4-rms, normalized) 

0.59 TT-mm-mrad 

0.39 TT-mm-mrad 


- 33 - 









0.08 



z (m) 

Figure 24. Horizontal beam envelopes of configuration I: 230 A on SNl. Protons are focused at the location 
of the 4-Jaw slits. Beam line elements are: ion source (VIS), solenoid 1 (SNl), segmented collimator (4Set), 
plain collimator (Cl), small-ID adapter tube (AT), movable slits (4Jaw), emittance scanner (ES), Faraday 
cup (FC). 


4.2 Comparison of WARP LEBT simulations with BCS results 
4.2.1 Simulations of Configuration I 

The first part of the simulation study was Configuration I, which included the beam line up 
to Position I, where a Faraday cup was located. This was also the first experimental configuration 
and direct comparisons of beam currents, emittance scans, and beam cross-sections are shown in 
this Section. The solenoid 1 (SNl) power supply was limited to a current of 350 A and the beam 
energy (ion source high voltage) was reduced to 55 keV in most experiments to be able to bring the 
HJ beam to a focus further upstream. Hence the simulations were performed at the same reduced 
energy. The first experimental study was to increase the current of SN1 to focus first the protons 
and then the Hj into the Faraday cup. The result of this scan, together with the simulated values, 
was shown in Figure The ratio of protons and in the simulation was obtained from the 
measured spectrum. Usually, the drain current on the source high voltage platform power supply is 
a good estimate of the total extracted current. In the simulation, the measured total extracted current 
had to be reduced from 29.6 mA to 20 mA in order to get good agreement. The reason for this is 
unclear at this point. Possible explanations are a erroneous reading of the drain current (caused by 
secondary electrons), a larger than simulated beam halo, scraped at the 4-sector collimator (4Set), 
losses through residual gas interaction, or a combination of underestimated divergence of the initial 
beam and overestimated space charge compensation in the simulation. Running the simulation 
with the reduced beam current yielded good agreement with the solenoid scan and also with the 
emittance scans, so it stands to reason that the results of this section are a good starting point for 
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Figure 25. Horizontal beam envelopes of Configuration I: 320 A on SN1. is focused at the location of 
the 4-Jaw slits. A strong proton focus at 1.2 m can be seen. This leads to the formation of the hollow beam. 
Beam line elements are listed in Figure 
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the subsequent simulations presented in the next section. The beam envelopes corresponding to 
the proton peak and the peak in Figure are shown in Figures ^ and ^ respectively. 
Approximately 2.6 m after the source are the 4-jaw slits (4Jaw) which were completely closed at 
one point, to obtain photographs of the beam cross-section (cf. Figure |^. A comparison can be 
made with the simulations by creating x-y density plots from the particle distributions at the same 
location in the two simulation sets for SNl = 230 A and SNl = 320A. This is shown in Figure 
PSl . Both images agree well with the photographs in both size and shape of the distributions (note 
the ring shape of the beam in the right plot), except for one detail: the halo around the 
proton beam in the left photo of Figure im already forms a ring even before the protons come 


Table 6. Simulated diameters and emittances for protons and 3.5 m downstream of the ion source 
(location of the Faraday cup). Values were obtained at the respective minimum diameters (SNl at 220 A for 
protons and at 340 A for Hj. 



protons 

H2+ 

Energy 

55keV 

55keV 

x-diameter (2-rms) 

21.2 mm 

44.1 mm 

y-diameter (2-rms) 

21.2 mm 

40.1 mm 

x-x’-emittance (4-rms, normalized) 

0.63 ;r-mm-mrad 

1.24 TT-mm-mrad 

y-y’-emittance (4-rms, normalized) 

0.66 ;r-mm-mrad 

1.14 TT-mm-mrad 
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Figure 26. Colored current density plots of the beam cross-sections at the location of the 4-jaw slits (compare 
to Figure [^. Left; 240 A on SNl, Right; 320 A on SNl. The ring structure of the hollow beam is 
clearly visible on the right. 


to a tight focus upstream of the 4-Jaw slits. This is a hint that space charge compensation could 
be overestimated in the simulation. The diameters and emittances of the proton and focus 
in the Faraday cup (FC) are listed in Table ^ While the proton emittances increase by 17% and 
12% (x-x’ and y-y’, respectively), the emittances increase by 244% and 192% (x-x’ and y-y’, 
respectively). This is due to the effect that the protons have on the ions. The severe over¬ 
focusing (seen in Figure hollows out the beam and leads to a tripling of the emittance. This has 
been seen in the measurements as well. A comparison of the horizontal phase spaces (measured 
and from the simulation) is shown in figure 0. Both plots exhibit similar size and structure. The 
proton intensity in the simulation seems to be higher and the asymmetry of the beam shifted 
left rather than right, which we attribute to steering effects. Overall, the agreement is good. 
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Figure 27. Comparison of measured horizontal phase space (left) and current density plot of simulated 
horizontal phase space (right) with SNl at 320 A. The ’N’-shape is typical for a hollow beam and leads to 
signihcantly increased emittance. The beam energy is 55 keV in this case. The different intensity scales 
arise from the fact that the entrance slit of the emittance scanner was smaller than the step size. 
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Figure 28. Horizontal beam envelopes of configuration II: 353 A on SNl, 230 A on SN2. is focused at 
the location of the Faraday cup at the end of the beam line which roughly coincides with the entrance aperture 
of the spiral inflector. Beam line elements are: ion source (VIS), SNl, segmented collimator (4Set), plain 
collimator (Cl), small-ID adapter tube (AT), movable slits (4Jaw), pneumatic beam stop + fixed collimator 
(Ch-BS), SN2, emittance scanner (ES), Faraday cup (FC). 

4.2.2 Simulations of Configuration II 

Here we present the results of the simulations of the full LEBT guiding the beam to the 
cyclotron. The particle distributions at the end of these simulations were used as initial conditions 
for the spiral inflector simulations reported in Section P3| . A typical horizontal beam envelope is 
shown in Figure As in Configuration I, the tight focal point of protons early on can be seen, 
which leads to sttong aberrations in the beam. Simulated phase spaces at the location of the 
emittance scanner ("‘ES"’ in Figure are shown in Figure ^ which should be compared to the 
measured phase spaces of Figure E Good qualitative agreement can be seen. Protons are faint 
streaks in the back, while the can be focused to below 11 mm 2-rms diameter. The diameters 
and emittances of the beam for different settings of SN2 are compared to the measured values 
in Figure [T^ It should be noted that, in the measured phase spaces, part of the beam is outside 
of the axes limits and thus the calculated values are underestimated. In Figures 0 and ^ the 
same limits were applied to the simulation results leading to good agreement. When using the 
full untruncated particle distributions to calculate diameters and emittances from the simulations, 
the diameters remain largely unchanged, while the emittances increase. The average increase of 
emittances from truncated to untruncated phase space is « 17.3% for the horizontal phase spaces 
and Ri 15.4% for the vertical ones, with a maximum of 35.6% and 31.4% for SN2 at 200 A (largest 
beam). This is treated as a systematic error and is included in the errorbars in Figure [T^ 
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Figure 29. Current density plots of simulated phase spaces (compare to Figure 
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Figure 30. The final particle distributions obtained by the WARP LEBT simulations, Left; Final beam 
cross-section. Right: Final x-x’ phase space. Red: Hj, blue: protons. 


Using the LEBT settings (solenoid currents and 4-Jaw slit openings) used during the cyclotron 
injection tests in the Configuration II simulation, we obtain realistic particle distributions for further 
simulation studies of the spiral inflector transmission and RF capture in the first three turns of the 
cyclotron. The cross-section and horizontal phase space are shown in Figure ^ and the important 
beam parameters are listed in Table 0. Results of OPERA ray-tracing calculations using these 
distributions are reported in the following section. 


Table 7. The final beam parameters for at the cyclotron entrance aperture for the beam injected during 
the tests at BCS in Vancouver. 



Energy 

62.7 keV 

x-diameter (2-rms) 

10.6 mm 

y-diameter (2-rms) 

10.2 mm 

x-x’-emittance (4-rms, normalized) 

1.19 TT-mm-mrad 

y-y’-emittance (4-rms, normalized) 

1.16 ;r-mm-mrad 

4-rms includes 

85% 
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Figure 31. VectorFields OPERA ray-tracing results with different dee voltages. Left: V(jee = 70 kV (nominal 
value). Right: Vdee = 45 kV (minimum simulated). It becomes immediately clear that even perfectly in 
phase, the low dee voltage is not sufficient to capture the particles and most of them are lost in the first turn. 
Coloring of the trajectories is loosely based on time-of-flight, but should be considered as a means to guide 
the eye only. 


4.3 Comparison of injection measurements with OPERA simulations 

As mentioned in Section [Z.1.3| , VectorFields OPERA was used for the design of the spiral 
inflector and the central region of the BCS test cyclotron. OPERA uses a finite elements method to 
calculate electrostatic and magnetostatic fields and has a ray-fracing module fo frack ions fhrough 
fhe simulafion space. In fhis section, we presenf a sef of ray-fracing calculafions of a 62.7 keV 
beam enfering fhe cyclofron fhrough fhe spiral infleclor and being accelerafed for 3.5 furns. 
The simulafion paramefers are lisfed in Table and are fhe same as used during fhe measuremenfs 
reported in Secfion ^31 The initial parficle disfribulion for fhe OPERA simulafion was obfained 
from fhe WARP EEBT simulafion discussed in fhe previous secfion. A unbunched beam was 
assumed and a phase resolufion of 5 degrees wifh 500 macro-parficles per phase was used during 
fhe OPERA calculafion. A variefy of 6 differenf dee volfages Vdee was used in order fo invesfigafe 


Table 8. Parameters for inflection and acceleration simulations using OPERA. 


Parameter 

Value 

Species 

H+ 

Initial Beam Energy 

62.7 keV 

Solenoids 1/2 

350 A/233 A 

Cyclotron Magnet 

218.6 A 

Spiral Inflector Upper / Eower Electrode 

-10.0kV/-tl0.15 kV 

Dee voltage 

45 - 70 kV 
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Figure 32. VectorFields OPERA calculation of transmission through the central region of the cyclotron. 
The counting of completed turns starts at the exit of the spiral inflector. As comparison, the measured 
transmission for reduced dee voltages of are also plotted. A clear reduction in transmission can be seen for 
lower dee voltages. 

the effect of insufficient dee voltage on the transmission (as encountered during the experiment). 
Figure ^ shows one sample of trajectories each for the cases of Vdee = 45 kV and Vdee = 70 kV, 
perfectly in phase with the cyclotron RF. 

The total transmission through the first turns of the cyclotron for different voltages applied to 
the dees can be seen in Figure ^ as well as the two measured curves. The initial current of the 
simulations was scaled to 5 mA which is the value assumed at the exit of the spiral inflector during 
the measurement. It is not surprising, that the transmission decreases for dee voltages below the 
design value of 70 kV. Vdee = 50 kV agrees reasonably well with the better of the two measured 
curves. The average dee voltages during the measurements determined by the pick-up probe and 
from beam radii (radial probe positions) were between 50 and 60 kV. Based on the large variation 
of the pick-up probe signal during the measurements, it is reasonable to assume that the average dee 
voltage of 48 — 50 kV giving good agreement with simulations is what could be achieved with the 
given experimental setup. Any discrepancy would be due to the high uncertainties in the pick-up 
probe values and the unstable dee voltage. The calculation of beam energy from the mean beam 
radii (radial probe positions) and magnetic field data (cf. Table ||) suggests similar values. 

Space charge effects, which could lead to additional reduction in transmission during the first 
turn, were not included in the OPERA calculations. This will be investigated in the future by 
using the particle-in-cell (PIC) code OPAL [35], which has a dedicated cyclotron module that has 
recently been upgraded to include spiral inflectors. 
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5. Conclusion 


This technical report outlined the experimental results obtained at the Best Cyclotron Sys¬ 
tem’s test bench, exploring design concepts and hardware relevant to the IsoDAR injection sys¬ 
tem. In addition we have benchmarked the results against a number of detailed simulation studies. 
These experimental results represent the early stages of exploration into the development of the 
IsoDAR injection system and provided excellent information regarding working with high-current 
HJ beams. Some specific highlights of our measurements, described in this paper are: 

1. The VIS source was able to produce good currents of H 2 ions (9.0 ±0.5 mA); with the initial 
(large) plasma chamber the best ratio of H^/protons was > 40%, obtained at microwave 
power just above the threshold for stable discharge operation (about 300 watts). This current 
improved to 12.2 ± 0.6 mA using a smaller (2.5 cm dis x 10 cm long) plasma chamber, 
and at the same time the H^/proton ratio improved to > 50%. Modeling and systematic 
measurements had suggested that smaller chambers would improve production, due to 
wall-recombination effects, and this was indeed borne out by experiment. We expect further 
optimization of the chamber shape to yield yet improved performance. 

2. Separation and elimination of the proton fraction from the H 2 beam was attempted by using 
focusing solenoids instead of magnetic analysis with a bending dipole, because of constraints 
in the existing test bench at BCS. This separation was successful in that very few protons 
were seen at the location of the test cyclotron, however as the separation process involved 
over-focusing the proton beam, which was then lost on the beam pipe walls, substantial 
heating of these walls was observed. In addition, the strong electric fields af fhe fighf proton 
focus had significanl disrupfive effecls on fhe H 2 beam. In parficular, our sfudies showed 
fhaf downsfream of fhis profon focus fhe Hj beam was hollow: ifs cenfral region had been 
blown ouf by fhe sfrong fields af fhe profon focus. We conclude fhaf fhe use of a dipole 
magnef would mifigafe such effecls and improve overall fransmission efficiency and 
beam qualify. 

3. The imporfance of space-charge compensalion in fhe Iransporl line was clearly observed, 
wifh subslanfial increases in beam size visually seen when elecfric fields swepf away fhe 
compensating elections. Measuremenfs wifh a relarding field analyzer were performed, and 
will be reporled in a subsequenl paper. 

4. A well-characlerized 6.2 ±0.4 mA, wifh measured normalized 4-rms emillance of 1.15 ± 
0.09 TT-mm-mrad af fhe enlrance lo fhe spiral infleclor was tiansmifled Ihrough fhe inflecfor 
wifh heller lhan 90% efficiency. This validales bolh fhe design and fabricalion of fhe spiral 
inflecfor. 

5. Beam was caplured and accelerafed Ihrough 4 lurns in fhe cyclotion, however wifh very high 
losses because of inadequate RF accelerating vollage. The problem of sparking across a 
critical insulalor which limiled fhe highesl vollage oblainable, and unforlunalely fhis insula¬ 
tor could nol be redesigned before fhe end of our experimenfal run. However VecforFields 
OPERA simulalions did provide qualifalive fils to fhe dafa observed. 
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6. A very successful set of simulation studies was able to reproduce quite well the observed 
measurements, including the pronounced emittance aberrations and hollow beams; indicating 
that the beam dynamics in our transport line can be understood and modeled with these 
well-established simulation codes. These then are now well validated for the design of the 
injection line for our future Hj experiments. 

Building on the base of these measurements, we see several possible clear paths forward for 
designing the IsoDAR injection line: 

• Ion source: The current generation of non-resonant microwave ECR sources, typified by 
the VIS we used, is probably capable of up to about 15 to 20 mA of Hj, when fully op¬ 
timized. This may be adequate for IsoDAR provided an efficient RF bunching system can 
be developed. While it was our initial aim to test such a bunching system, the problems 
with the cyclotron RF system prevented the installation and running of the buncher that had 
been built. If a bunching efficiency of 20 fo 25% can be achieved, the 5 mA injected current 
requirement could be met. We will endeavor to find a way of conducting these tests in the 
future. 

• Inflection: A key result of our experiments has been the validation of the spiral inflector 
design for these more rigid, high current beams. It should be noted here that the final IsoDAR 
design will be a further improvement over the BCS test stand model as it will operate in 
combination with 4 RF cavities in 4* harmonic mode. This will allow for better vertical 
focusing and will improve the injection efficiency of the cyclotron. 

• Alternate design: A new FEBT concept incorporating a low-frequency RFQ (Radio Fre¬ 
quency Quadrupole) structure as a buncher might be an elegant, compact and efficient means 
of providing beam to the IsoDAR cyclotron. Such a device [36], operating at the cyclotron 
frequency and placed in the axial bore hole just above the center of the cyclotron could pro¬ 
vide in principle bunching of up to 80 or 90%, greatly relaxing the requirements on source 
performance. In addition, the injection line could be much reduced in size and length, which 
might also potentially lead to overall cost savings. Initial studies are underway [37]. 
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